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Abstract. We discuss multi-hadron production in both inelastic neutrino-nucleon interactions in the cur- 
rent fragmentation region and neutrino-nucleus collisions in the target fragmentation region. Our analysis, 
carried out within the framework of the quark-gluon string model, is mainly focused on the difference 
between these two processes. We show that the Q 2 dependence of hadron multiplicity in the current and 
target fragmentation regions is indeed completely different. The study of inelastic v — A scattering in the 
target fragmentation region also provides new information on nuclear structure at small N — N distances. 
The results of the proposed approach are in satisfactory agreement with the data recently obtained at 
CERN by the NOMAD Collaboration. 

PACS. 13.60.Hb Total and inclusive cross sections - 13.15.+g Neutrino interactions 
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1. Introduction 

Deep-inelastic particle interactions play a decisive role in 
the development of the modern theory of fundamental in- 
teractions. For instance, the investigation of multi-hadron 
production in inelastic scattering of leptons off nucleons 
and nuclei is a tool to study both the dynamics of such 
processes and the quark structure of the target. Of par- 
ticular interest, in this context, is the analysis of the Q 2 
dependence over a wide range of values. The difficulties as- 
sociated with the increase of the QCD coupling constant 
a s (Q 2 ) at moderate and low Q 2 have led to the devel- 
opment of alternative approaches to analyze soft hadron 
interactions. A well known example is the the 1/N ex- 
pansion of the scattering amplitude, N being the number 
of colors or quark flavors, suggested by t'Hooft P and 
Veneziano Based on this approach, the Quark-Gluon 
String Model (QGSM) and the Dual Parton Model (DPM) 
have been developed in Refs.jS] and @E], respectively, to 
analyze hadronic reactions. The first application of the 
QGSM to inelastic lepton-nucleon processes at moderate 
and low Q 2 has been carried out in Ref.|S]. 

In this paper we investigate inelastic neutrino scatter- 
ing from both protons and nuclei within the QGSM, fo- 
cusing mainly on the differences between these reactions. 
In Section 2 we outline the main features of the theoreti- 
cal approach, whereas Section 3 is devoted to a compari- 
son between our results and the data recently obtained at 
CERN by the NOMAD collaboration. Finally, in Section 
4 we summarize our findings and state the conclusions. 



2. Theoretical framework 

In the QGSM pOE], hadron production in the reactions 
v{v) + p — > + h + X is described in terms of pla- 

nar and cylindrical graphs, as shown in Fig^ The planar 
graph of panel (a) describes neutrino scattering off a va- 
lence quark, corresponding to one Reggeon exchange in 
the ^-channel jSJ, whereas the cylindrical graph of panel 
(b) describes neutrino scattering off sea quarks, corre- 
sponding to onc-Pomeron exchange in the t-channel 0. 
The figure also shows the occurrence of hadronization in 
the colorless quark-antiquark and quark-diquark strings. 

The relativistic invariant distribution of hadrons pro- 
duced in the process v(v)+p — > + h + X is defined 
as 
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where Eh and p h are the total energy and three momen- 
tum of the produced hadron, respectively, whereas E' and 
fi are the energy and the solid angle of the final state 
muon. The right hand side of Eq. (TjJ can be written in the 
following general form |7||S]: 
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Fig. 1. Planar, one-Reggeon exchange (panel (a)), and cylindrical, one-Pomeron exchange (panel (b)), graphs (see text). 



where G is the Fermi weak coupling constant, E is the en- 
ergy of the incoming neutrino, m and mw are the nucleon 
and PF-boson masses, respectively, x = Q 2 /2(p l/ ■ k) is the 
Bjorken variable and p v and k are the four momenta of 
the initial neutrino and nucleon. In general, z is the light 
cone variable defined as z = (Eh + ptz) I (E + p z ), where 
Phz is the component of the hadron momentum parallel 
to the momentum of the incoming neutrino, p z , (see, e.g., 
Refs.|HK] and 0JE1)- At large energies of the final hadron 
it reduces to the longitudinal momentum fraction of the 
hadron, with respect to the neutrino, in the rest frame 
of the target proton. The variable z can also be treated 
as the Feynman variable xf = 2p* L /Wx, defined as the 
longitudinal momentum fraction in the hadronic center of 
mass system (HCMS). Here p* L denotes the longitudinal 
hadron momentum in the HCMS and Wx is the mass 
of the hadrons produced in the reaction. Finally, pht is 
the transverse momentum of the produced hadron with 
respect to the current (hadronic jet) direction, (see, e.g. 
Ref.0). 

The probabilty distributions of hadron production as- 
sociated with the planar and cylindrical graphs of Fig^ 
are given by F P (x,Q 2 ; z,p ht ) and F c (x,Q 2 ; z,p ht ), re- 
spectively. They can be computed cither analytically, as in 
Ref.[7j, or using the Monte Carlo (MC) approach to gener- 
ate all quark-antiquark and quark-diquark strings shown 
in Fig^ as done in Ref.jH]. 

The main ingredients for the calculations of observ- 
ables in the reaction under discussion are the quark dis- 
tributions in a nucleon and their fragmentation functions 
to hadrons. In addition to the dependence upon x, quark 
distributions also depend on Q 2 and the transverse mo- 
mentum k t . Following Ref. J2j, we use a factorized form 
for these distributions: 



q f (x,Q 2 ;k t ) = qf(x,Q 2 )g q (k t ) 
with the function g q chosen in the form 
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where B = l/(k t ) ~ 4(GeF/c)- 1 , while (k t ) ~ 0.25 GeV/c 
is the average transverse momentum of a quark in a nu- 
cleon. As for the function qf(x,Q 2 ), we use the fit sug- 
gested in Ref. 10 , including true Regge x-asymptotic at 
x — > 0, x — > 1 and small Q 2 , and its QCD prediction at 
large Q 2 . 

In general, the fragmentation functions (FF) of quarks 
(diquarks) into hadrons, , depend on the hadron mo- 

mentum fraction z\ and the hadron transverse momentum 
with respect to a quark (diquark) momentum direction, 
Pht. Here we choose again the factorized form 



D q(qq)( Z l,Pht) = D^ qq) (Zl)g q (p ht ) , 



(6) 



with the function g q defined as in Eq.JSJ- 

The functions D q ,. (z%) are obtained, according to the 
recursive cascade model procedure suggested in ^T] , from 
the integral equation 
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the function f(x) being chosen in the form 

fix) = x P(i-xy . (8) 

According to the main assumptions underlying the QGSM, 
the fragmentation functions D q ( qq )i z i) should satisfy true 
Regge asymptotic at z% — ► 1 and z\ — > 3.. These con- 
straints determine the values of the parameters. The de- 
tailed procedure is presented in Ref.^2]- I n general the 
FF depend not only on z\ and pht, but also on Q 2 . At low 
Q 2 they have to reproduce the true Regge asymptotic |3] , 
while at large Q 2 they have to describe e + e~ annihilation 
data. As we mainly analyze inelastic v — N interactions at 
moderate Q 2 , one can assume that the Q 2 -dependence of 
the FF be negligibly weak. 

Let us now consider pion production in the target frag- 
mentation region in v — A interactions. Our study of these 
processes focuses on the kinematical region of large pion 
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momentum (p„ > 0.3 GeV/c), where the effects of final 
state interactions (FSI), leading to pion absorption asso- 
ciated with production of baryon resonances, are expected 
to be small |13U14| and the impulse approximation (IA) 
cane be safely used. 

Within the framework of the IA the relativistic in- 
variant semi-inclusive spectrum of hadrons, in particu- 
lar pions, produced by the semi-inclusive inelastic process 

e + A^'e + tt + x, 



Pe+A^e'+TT+x = Ehd- 



(PpudQdE 1 
can be written in the convolution form [7] 



(9) 



e+A^e'+7T+x(x,Q ;z,p t ) = / dyd k t fA(y,Q ,h) 

Jz<y 
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where Z , N and A are the numbers of nucleons, protons 
and neutrons in the nucleus, z = (pp u ) / 1 (PaPv)Mai 'm, 
p is the four-momentum of the produced pion and Ma 
and m arc the nucleus and nucleon masses, respectively. 
Pt+p(n)— >t'+ir+x denotes the semi-inclusive spectra of pi- 
ons produced by interactions of the lepton £ with a quasi- 
free proton (neutron). The nucleon distribution function 
kt) is defined as JS] 



./ ,(//■ h I = / dk dk z S(k)yS (y - ^-^y 



, (11) 



where q is the four-momentum transferred by the lepton, 
q 2 = — Q 2 , S(k) is the relativistic invariant function de- 
scribing the nuclear vertex with an outgoing virtual nu- 
cleon, y = (Mji/ni)(kpi,) / (PaPv) and Pa >k andp„ are the 
four-momenta of the nucleus, nucleon and initial neutrino, 
respectively. 

The distribution function /^(y), defined as the integral 
of fA(y,kt) over d 2 k t , can be calculated within nuclear 
many-body theory approximating S(k) with the nonrela- 
tivistic spectral function P(fc, E), yielding the probability 
of finding a nucleon with momentum k and removal en- 
ergy E — m—ko 46 . However, due to the limited range of 
momentum and removal energy covered by nonrelativistic 
calculations of P(k,E) (typically |k| < k m ; m ~ 0.7 — 0.8 
GeV/c and (m — k ) < 0.6 GeV, see e.g. JE]), this proce- 
dure can only be used in the region y < y^ ~ 1.7 — 1.85. 
An alternative approach to obtain /a(?/) at larger y, based 
on the calculation of the overlap of the relativistic invari- 
ant phase-space available to quarks belonging to strongly 
correlated nucleons, has been proposed in Refs. |15U17| . 

The convolution form of Eq. (|ll|l implies that to ana- 
lyze the production of pions emitted in the whole back- 
ward hemisphere by v — A interaction one needs to know 
the dependence of the nucleon distribution upon y and kt ■ 
We assume the following factorized form: 



where the function qa is chosen to be a Gaussion 
9A(k t ) = — T2rexp(-fc t 2 /(fc t 2 )) , 



(13) 



(fc 2 ) being the average value of the squared transverse mo- 
mentum of a nucleon in the target nucleus. 



3. Comparison to data 

The results of our calculation of different observables of 
the reactions v + N^fi~+h + X and v + C 12 — > fi~~ + 
7r~ + X are presented in FigsGhE] 

The mean charged multiplicity in the current region 
is shown in Fig|3 The multiplicity measured by the NO- 
MAD collaboration [T^], (n c h), turns out to be close to 
(n c /j)/2 measured in e + e~ experiment at energy E — ^/s 
and to (n c h) obtained from ep and v^p at E — Q. In Fig|2| 
we compare our results to (n7 h D )/2 obtained from a QCD 
calculation of the charged multiplicity in e + e~. The QCD 
results, represented by open circles, correspond to evolu- 
tion of parton distributions in the leading log approxima- 
tion, yielding the following fit for (n^ D ) (see Ref.fTH] 1 ): 



QCD 



b exp cJlog(Q 2 /Q 



(14) 



where a = 2.257, b = 0.094, c = 1.775 and Q = 1 GeV/c. 
The solid, long-dash and short-dash lines in Fig[5] corre- 
spond to our Monte Carlo calculations at initial neutrino 
energies E v = 150 GeV, E v = 45 GeV, E v = 23.6 GeV, re- 
spectively. The NOMAD experimental data are averaged 
over the initial energy |18j. 

FigEl shows that the QCD fit, while failing to repro- 
duce the NOMAD data on the charged multiplicity at 
Q < 5 GeV/c, provides a reasonable description of the 
behavior at large Q. As discussed in ^Sj the same fit also 
satisfactorily describes (n c h) obtained from e + e~ , ep and 
v^p experiments. On the other hand, the approach pro- 
posed in Ref.[5] describes the NOMAD data at Q < 5 
GeV/c rather well. 

The multiplicity of pions, normalized to the cross sec- 
tion of the process v + A — » /i + X, a, is defined as 



1 

a 

dz 

z 



dx 

d 2 p Pv+A^h+ti+X 



(15) 



fA(y,h) = fA{y)9A(k t ) 



(12) 



It can be split in two parts, corresponding to the planar ,or 
one-Reggeon exchange, diagram (Fig^(a)) and cylindri- 
cal, or one-Pomeron exchange graph (Fig^(b)). 

The NOMAD collaboration carried out a study of in- 
elastic v~C interactions in which the negative pions emit- 
ted in the whole backward hemisphere, with respect to the 
incoming neutrino beam, were detected The data of 
FigEl show the nultiplicity of negative pions carrying mo- 
menta 0.35 < p,r < 0.8 GeV/c, measured in a kincmatical 
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Fig. 2. Mean multiplicity of charged hadrons in the current 
fragmentation region for the process v + p — > [>T + h + X, 
as a function of the momentum transfer Q. The open circles 
correspond to the QCD fit given by ea. l|14[l . while the solid, 
long-dash and short-dash lines correspond to our calculations 
at E v = 150, 45 and 23.6 GeV, respectively. The experimental 
data are taken from Ref. |18| 
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Fig. 3. Mean multiplicity of charged pions produced in the 
backward hemisphere in v — C interaction as a function of Q 2 . 
The solid line shows the results of the full calculation, includ- 
ing both graphs of Fig^ The dashed and dash-dot lines corre- 
spond to the separated contributions of the cylindrical graph 
of FigQ(b) and the planar graph of FigQ(a), respectively. The 
dots show the results obtained setting P(k, E) = in the high 
energy-momentum domain, not covered by the nonrelativistic 
calculation of Ref.|lri). The experimental data are taken from 
Ref.0. 
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Fig. 4. -dependence of the spectrum of backward pions 
produced in inelastic v — C interaction. The solid curve corre- 
sponds to the full calculation whereas the dashed line has been 
obtained including only the cylindrical graph of FigQ(b). The 
dash-dot line shows the results obtained setting P(k, E) = in 
the high energy-momentum domain, not covered by the non- 
relativistic calculation of Ref.|lri). The experimental data are 
taken from Ref.pj- 



setup in which Wx increases as Q 2 increases. Theoreti- 
cal calculations have been performed applying the same 
kinematical conditions. 

As shown in Fig|3 the planar graph provides the main 
contribution at small Q 2 , while the cylindrical one domi- 
nates at Q 2 > 10 (GeV/c) 2 . It clearly appears that both 
diagrams have to be included to explain the observed Q 2 
dependence. The dotted line in FigEl has been obtained 
setting P(k, E) = in the domain of large energy and 
large momentum, not covered by the nonrelativistic calcu- 
lations of ref. ^1 ■ Comparison between the solid and dot- 
ted line shows that the dominant contribution to < > 
does indeed come from the high momentum tail of the 
nuclcon distribution. 

In FigQ]we compare the p 2 dependence of the experi- 
mental spectrum to the one resulting from our approach. 
The contribution of the cylindrical graph dominates the 
spectrum and the inclusion of the high momentum tail of 
the nucleon distribution, corresponding to p > 0.4 GeV/c, 
is clearly needed to describe the data. 



4. Summary and conclusions 

In this paper we have analyzed inelastic v — N and v — A 
processes at low and moderate Q 2 within the QGSM. The 
main conclusions can be summarized as follows. The dy- 
namics of hadron production in the current fragmentation 
and target fragmentation regions is different. At low Q 2 
the contribution of the cylindrical graph to the hadron 
multiplicity is much more prominent in the first kinemati- 
cal region than it is in the second one. There is some anal- 
ogy between pion production in the backward hemisphere 
from target fragmentation and soft h ~ N processes at 
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xf —* — 1, where the one-Pomeron exchange graph domi- 
nates at large initial energy (see, e.g., Ref.[H]). In conclu- 
sion, the contributions of the graphs of FigsQ]( a ) an d (b) 
have different Q 2 dependences in the different kinemati- 
cal regions corresponding to current and target fragmen- 
tation. 

The conventional perturbative QCD calculations of ha- 
dron multiplicities, both in v — N and v — A inelastic in- 
teractions, do not reproduce the NOMAD data, whereas 
the application of the QGSM allows for a rather satisfac- 
tory description. As for the role of nuclear structure, the 
dominant contribution to the spectrum of pions emitted 
in the backward hemisphere in v — A inelastic scatter- 
ing is clearly coming from the high momentum tail of the 
nucleon distribution, which can be described in terms of 
overlaps of distributions of three-quark colorless objects 
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